Evidence for Disorder Induced Derealization in Graphite 
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We present electrical transport measurements in natural graphite and highly ordered pyrolytic 
graphite (HOPG), comparing macroscopic samples with exfoliated, nanofabricated specimen of 
nanometer thickness. Nano-samples exhibit a very large c-axis resistivity p c - much larger than 
expected from simple band theory - and non-monotonic temperature dependence, similar to macro- 
scopic HOPG, but in stark contrast to macroscopic natural graphite. A recent model of disorder 
induced derealization is consistent with the data. Micro-Raman spectroscopy reveals clearly re- 
duced disorder in exfoliated samples and HOPG, as expected within the model - therefore presenting 
further evidence for a novel paradigm of electronic transport in graphite. 



Graphite is a paradigmatic layered material and has 
been investigated intensively for many decades. The in- 
plane resistivity p a b is rather well described by a simple 
Drude model. However, the resistivity p c along the c- 
axis, perpendicular to the graphite basal planes, as well 
as its temperature dependence p c (T), are not described 
by the simple band structure model [TJ [2] , and currently 
lack theoretical understanding despite extended efforts. 
The resistive anisotropy Ra = p c /Pab is a convenient 
dimensionless parameter characterizing transport prop- 
erties. 

Carbon atoms in the graphite basal planes are strongly 
bound by covalent bonds, while much weaker Van der 
Waals forces bind the graphene sheets along the c-axis. 
Non-trivial disorder such as stacking faults and crys- 
talline grains result in a mosaic angle and complicate 
electronic transport. Simple band theory [F predicts 
Ra — m c /m a 5 ~ 140, the ratio of the corresponding 
band masses, which holds in presence of isotropic dis- 
order. This agrees well with measurements in natu- 
ral graphite (NG) [3j 0]. In highly oriented pyrolytic 
graphite (HOPG), the anisotropy was found to be much 
larger, even exceeding lO'OOO in some experiments [5j|6]. 
Moreover, band theory pQ predicts a monotonic metal- 
lic temperature dependence for both p^ and p c , result- 
ing in a temperature independent anisotropy Ra- This 
is seen in NG [7], but not in HOPG, where p c is non- 
monotonic with a maximum around 40 K [5, 6, 8 , simi- 
lar to p c in other layered materials, such as NaCo2 04 [9] 
and Cuprates [10]. A large anisotropy far exceeding 100 
combined with the non-metallic temperature dependence 
- together referred to as anomalous behavior - are cur- 
rently not understood and present a fundamental prob- 
lem in condensed matter physics [7J [TTJ [12] . 

In this Letter, we use exfoliation and nano-fabrication 
techniques to investigate electrical transport in graphite 
flakes (see Fig.JT]), comparing samples with thicknesses 
in the nanometer range with macroscopic samples and 
further comparing HOPG with Indian and Madagascar 
NG. Surprisingly, while conforming to the band-structure 
results in macroscopic NG flakes, NG displays anoma- 
lous behavior in nanostep samples - namely a large re- 



sistive anisotropy as well as a non-metallic temperature 
dependence. HOPG, on the other hand, shows no size 
dependence, behaving anomalously in all samples. We 
note that there are no previous p c data for nanometer 
thickness specimen, to the best of our knowledge. Inter- 
estingly, the in-plane resistivity of all samples is in good 
agreement with reported values [3H5] , and shows no size- 
dependence, making it clear that the large anisotropics in 
the anomalous samples are to be attributed to a large p c . 
These measured anisotropics appear consistent with a re- 
cent model based on disorder induced derealization by 
Maslov et al. [13] [14] , further corroborated by a disorder 
characterization of our samples using micro- Raman spec- 
troscopy. Conduction path mixing due to a finite mosaic 
angle can account for the non-monotonic temperature de- 
pendence [15J, altogether presenting first experimental 
evidence for a novel paradigm of electrical transport in 
graphite. 




FIG. 1: Nano-graphite samples, (a) Device schematic. Ti/Au 
contacts (yellow) for 4-wire measurements are patterned on 
each plateau, isolated from the graphite walls by a S1O2 layer 
(purple). AFM picture (b) and optical microscope image (c) 
of an HOPG flake with two plateaus, (d) Cross section along 
the blue line in (b), giving plateau heights. 
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To produce nanostep samples, we use the design shown 
in Fig.[TJa). We exfoliate graphite onto a SiC^/Si wafer 
and identify suitable flakes with two plateaus differing in 
height using an optical microscope. The lower plateau 
height d and the step height h are determined from an 
AFM image and cross-sections, see Fig. [TJb,d), giving 
heights between 14 and 150 nm. To extend the range 
to larger step heights, we use e-beam lithography and 
oxygen-plasma etching to carve steps up to h = 450 nm. 

For contacting the plateaus, we first cover parts of 
the exterior edges of both plateaus with Si02 of at least 
80 nm thickness [16] in order to prevent short-circuiting 
the c-axis. Contacts (typically a few hundred nanome- 
ters in width) and bonding pads are patterned in a fi- 
nal e-beam step, evaporating a Ti/Au layer thicker than 
110 nm (Si02 thickness plus 30 nm). A typical device is 
shown in Fig. [ljc). All resistances are measured with 
a standard 4- wire lock-in technique. This layout allows 
measurements of both p^ and p c on the same device, 
as needed to obtain the anisotropy, but requires resistiv- 
ity corrections due to geometry inherent conduction path 
mixing. 

The in-plane resistivity p a b = R a b^ab/lab is estimated 
from the 4-wire resistance R a b with current and volt- 
age probes on the same plateau and assuming a simple 
rectangular shape of the graphite sample, with voltage 
probe distance l a b and total graphite cross section 
(see Fig.JTJa)). This is a good approximation for thin 
samples, if the anisotropy is smaller than l a b/d ~ 250. 
For large anisotropy - as seen in most samples - the ex- 
tracted p^ presents an upper bound due to an effectively 
thinner conduction layer, since the current cannot pen- 
etrate easily along the highly resistive c-axis. The p a b 
extracted here (see Table|lJ and Fig.|2| agree rather well 
with literature H7]. typically giving p^ ~ lpVtm. 
Moreover, p^ appears independent of the graphite thick- 
ness and is similar for NG and HOPG samples, as seen 
in Fig.[2j open symbols. 

Next, we determine the c-axis resistivity p c . Since Z c , 
the contact to contact distance across the step, is much 
larger than the step height h (see Fig.JTJa)), we need to 
subtract the in-plane contributions to the measured re- 
sistance R c to obtain the corrected resistance R cor r using 
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with upper/lower contact to step distance l c i/ cu and cor- 
responding plateau widths wi/ u . p c is then obtained from 
p c = R corr A u /h, where A u is the upper plateau area. 
We note that as previously for p aoi we again overesti- 
mate the thickness d for large anisotropy. However, here, 
this effectively cancels the overestimated p aoi making the 
extracted p c quite robust. The correction is often sub- 
stantial, reducing p c and the anisotropy by a large factor 
(see Table]!]). The resulting p c are displayed in Fig. [2] as 
a function of height (filled symbols), giving very large 



p c and correspondingly large anisotropy Ra for all nano- 
graphites, both NG and HOPG. Interestingly, a power- 
law fit (linear fit on the log- log graph, slope — 1 ± 0.4) 
through all NG nanostep p c data points seems to indi- 
cate a trend reducing p c with increasing step height to- 
wards the macroscopic p c value in NG samples. HOPG 
nanostep data is excluded from the fit, since HOPG has 
no apparent size dependence when going from macro- 
scopic to nanostep samples (filled red diamonds). How- 
ever, the wide spread of the data and the resulting error 
bars on the fit-slope are significant; to make a more def- 
inite statement, samples with step heights between 1 pm 
and 100 pm might give more insight [T8] . 
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FIG. 2: Influence of graphite thickness on p c (solid markers) 
and pab (empty markers) at room temperature, comparing 
HOPG (red) with Madagascar NG (green) and Indian NG 
(blue). For p a b, the abscissa value used is d + h, the overall 
flake thickness, see Table|l] Previous measurements of macro- 
scopic samples (black) were added for both HOPG 5_ (stars) 
and NG (3, 4 (circles) for comparison. Dashed horizontal 
lines indicate literature values p ao ref for p a b and /? c ref for 
p c . Further, the best power-law fit to all NG nanostep data 
(dotted line, linear fit on log data, slope of — 1.0±0.4) is added 
to indicate a potential trend, see text. 

To allow a comparison with previous experiments, we 
also investigate macroscopic NG and HOPG samples, 
again measuring both p a ^ and p c . Due to the geome- 
try used, corrections are small and not necessary for the 
macroscopic samples. The values obtained are also added 
to Fig.[2j together with typical values from literature [3]- 
E] . We find decent agreement between our macroscopic 
data and previous measurements, reproducing here again 
the large discrepancy in p c between HOPG and NG in 
macroscopic samples. 

Next, we turn to the temperature dependence p c (T) 
of the macroscopic samples [19] . For HOPG, we find 
a non-metallic p c at high T (dp c /dT < 0), see Fig. [3] 
(a). Around 40 K, p c displays a rather shallow maximum, 
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TABLE I: Nano-sample parameters, see text for definitions, and Fig. [T] for an illustration. Samples listed here are represented 
in Figs.|2]and[3] Fig. [2] includes data from additional exfoliated Madagascar samples. 



Material 


p ah [pVtm] 




l c i[pm] wi[fim] I 


zu [pm] 


w u [pm] i 


i[nm] 


h[nm] 


R C orr[Q] A u [pm 2 ] p c [mQm] 


Ra 


HOPG 1 


1.22 


11.5 


4 


25 


2.5 


25 


30 


14 


2.3 


361 


58.0 


47700 


HOPG 2 


0.78 


12.1 


6 


20 


2 


20 


30 


24 


2.9 


57 


6.88 


8840 


Madagascar 1 


1.45 


27.4 


1.2 


10 


2.3 


10 


16 


80 


13.1 


248 


40.5 


27900 


Madagascar 2 


0.60 


26.0 


4.2 


8 


3.3 


5 


14 


37 


5.6 


95 


14.5 


24000 


India 1 


0.71 


18.9 


1 
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15 


17 


10.7 


100 


63.2 
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India 2 


5.27 


37.0 


2.4 


20 
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15 


20 


150 


1.27 


88 


0.78 


142 
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FIG. 3: Temperature dependence of resistivities. p c (T) in 
macroscopic HOPG (a) and macroscopic Indian (b) NG. (c) 
p a b(T) in nanoscale samples for HOPG (red) and NG India 
(blue) and Madagascar (green). Two samples are presented 
for each graphite type, corresponding to Tablejl] (solid, sample 
1; dashed, sample 2) (d) p c (T) for the same samples. 

in good agreement with previous HOPG measurements 
[5]. In contrast, macroscopic Indian NG behaves weakly 
metallic and monotonic down to 4 K (shown in Fig.|3](b)), 
also in agreement with previous NG data [3 j. Overall, 
our data from macroscopic samples fully agrees with the 
literature, giving us confidence that a comparison of the 
exfoliated samples with literature is appropriate. 

The temperature dependence of the exfoliated nano- 
graphite samples are shown in Fig. [3] (c) and (d), nor- 
malized to the high-T value. In most samples, p a b(T) 
is metallic and monotonous, as expected, and in agree- 
ment with macroscopic data |3j|20]. In two specimens, 
p a b is weakly non-monotonic. This seems to occur oc- 
casionally in nanoscale samples, as previously reported 
[2TI [22] • Interestingly, p c (T) of all nanostep samples is 
qualitatively the same, showing a non- metallic and non- 
monotonic temperature dependence, qualitatively iden- 
tical to macroscopic HOPG, and clearly different from 
the macroscopic NG data. We emphasize that the non- 
metallic p c (T) combined with the large anisotropy Ra 
constitutes anomalous behavior for all nanoscale samples. 
In contrast, only macroscopic HOPG is anomalous, not 



macroscopic NG. 

Motivated by an anisotropy far exceeding the band 
structure expectation, we consider a recent theory by 
Maslov et al. [13]. A similar effect was also previously 
observed for photons [13] , Within this theory, c-axis 
transport is strongly suppressed in samples with weak 
bulk disorder due to ID Anderson localization along the 
c-axis induced by randomly spaced barriers (e.g. stack- 
ing faults). This gives a very large p c and anisotropy 
Ra, as observed in the anomalous samples. In samples 
with sufficiently strong bulk disorder, however, c-axis lo- 
calization is destroyed by bulk scattering out of the c- 
axis direction, leading to reduced p c and smaller Ra, as 
seen in the macroscopic NG samples. Interestingly, here, 
disorder can suppress Anderson localization, rather than 
enhancing it, as is usually the case. Therefore, HOPG 
and nanostep samples are expected to have weak bulk 
disorder, whereas macroscopic NG specimen either have 
significantly more bulk disorder (suppressing c-axis local- 
ization), or fewer c-axis barriers, insufficient for localiza- 
tion (barriers spacing exceeding coherence length). 

To characterize disorder, we have measured spatially 
resolved micro-Raman spectra, see Fig. [4] We plot f = 
Id/Ig, the ratio of the G-peak intensity Iq and D-peak 
intensity Id after background subtraction, obtaining an 
intensity independent characterization of disorder. For 
graphite, < £ < 1. A large £ indicates a high degree of 
disorder (large Id). Due to a finite integration time, £ 
has measurement noise of ~ 1/50, i.e. £ < 1/50 indicates 
an invisible D-peak and weak disorder. We note that for 
macroscopic samples, the Raman measurements probe 
only a thin surface layer. For the nanostep samples, on 
the other hand, a significant depth is characterized [23J, 
presenting a good measurement of graphitic bulk disor- 
der. 

In macroscopic and nanostep HOPG, disorder is very 
l° w (£ ^ 1/50 essentially everywhere, not shown), while 
macroscopic Madagascar NG samples are clearly more 
dirty, consistent with expectations from the model. Mi- 
cron sized regions of strong disorder mix with cleaner 
areas in pristine NG, see Fig.[4^a) for a typical scan. A 
histogram is adjacent to the scan, revealing a broad dis- 
tribution extending up to £ ~ 1/3. In contrast, an ex- 
foliated Madagascar flake displays significantly reduced 
disorder, see Fig.Qb). Several exfoliated NG samples 
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FIG. 4: Disorder characterization with micro- Raman spec- 
troscopy. Ratio £ = Id/Ig of the D-peak and the G-peak 
intensities, scanned (resolution ~ 0.5 /xm) on (a) macroscopic 
Madagascar NG and (b) exfoliated Madagascar flake, plot- 
ted on the same color-scale, as indicated. White is off the 
graphite flake. Corresponding histograms of £ gathered over 
the shown scan areas are shown to the right of each scan. 

were inspected, all exhibiting £ significantly smaller than 
in macroscopic NG samples. 

We note that exfoliating macroscopic NG only once 
already results in strongly reduced disorder. Clearly, ex- 
foliating somehow acts expose clean surfaces [24] and to 
remove bulk disorder, though the mechanism is not clear. 
This could potentially explain the size dependence men- 
tioned in Fig.|2j thinner samples tend to require more 
exfoliation steps, therefore becoming cleaner, more lo- 
calized, and obtaining a larger p c . Ultimately, for suffi- 
ciently small /i, ID localization should break down and 
p c decrease strongly - not visible in the present data, 
presumably because h is still too large. Overall, the Ra- 
man data is in good agreement with the predictions of 
the model, namely weak disorder in all exfoliated and 
all HOPG samples, as opposed to strong disorder in NG 
macroscopic samples. 

Finally, we turn to the anomalous temperature depen- 
dence of p c . If the c-axis resistivity p c is very large, the 
c-axis conductance path could easily be mixed with the 
ab-conductivity a a b due to the finite mosaic angle ef- 
fectively short circuiting the intrinsic c-axis conductivity 
a c . Assuming small tilting 6 <C 1, the measured conduc- 
tivity <j c can be written as [15] 

a c (T) = a c (T) + (0 2 )-a ab (T), (2) 

where (6 2 ) is the variance of 6. In low bulk-disorder sam- 



ples at low temperatures, the intrinsic a c is very small 
(strongly localized) and a c (T) obtains a significant com- 
ponent from a a &, including the (weakly) metallic tem- 
perature dependence <r a 5(T), leading to a slight increase 
of p c (T) upon increasing T. At higher T, localization 
is weakened (due to phonon scattering, equivalent to in- 
creasing bulk disorder for increasing T), a c is enhanced 
and becomes increasingly more dominant, leading to a 
decreasing p c above some cross-over T. For a mosaic 
angle of about 0.8° (values between 0.2°and 2°were mea- 
sured, not shown), the crossover would occur around 
40 K, in agreement with the data. For disordered sam- 
ples, on the other hand, the intrinsic a c is dominating 
<J C (T) since localization is already lifted by disorder, re- 
sulting in the usual metallic temperature dependence, as 
seen in macroscopic NG (TJ- 

In conclusion, we observe anomalous behavior, namely 
high p c and non-monotonic p c (T), in both NG and 
HOPG exfoliated samples. This is in stark contrast to 
macroscopic samples, where the anomalous behavior is 
only seen in HOPG, consistent with previous experi- 
ments. A recently proposed transport theory [13] can 
consistently explain this convergence on the nanoscale, 
the macroscopic data, and the temperature dependence. 
Further, finding indeed low disorder in exfoliated and 
HOPG samples, and high disorder in macroscopic NG, 
as expected within the model. We note additionally that 
neutron irradiation experiments [5j [TT] inducing bulk dis- 
order also give consistent results, namely reduced p c after 
irradiation of HOPG, further corroborating the model. 
We therefore present first, clear evidence of disorder in- 
duced derealization as a new paradigm of electronic 
transport in graphite. 

Though beyond the scope of the present work, it 
would be very interesting to subject the model to further 
scrutiny: studying intermediate steps filling the thick- 
ness gap in Fig. [2] but also even smaller thicknesses, ul- 
timately down to few- or bi-layer graphene, potentially 
revealing the localization length. This might be facili- 
tated by bottom contacts with layers deposited on top, 
followed by top contacts. We note that the minimum 
thickness in the present samples is 14 nm, corresponding 
to about 50 graphene layers. Further, a characterization 
of graphite disorder would be of great interest, e.g. inves- 
tigating stacking faults and angles, intercalation, grain 
and boundary formation, aiming at identifying the local- 
ization mechanism, leading ultimately to a microscopic 
understanding of electrical transport in graphite. The 
results presented here were obtained in graphite, but it 
would be intriguing to learn if similar arguments apply 
to some of the numerous other layered materials. 
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